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The 21-cm anisotropies from the neutral hydrogen distribution prior to the era of reionization is a 
sensitive probe of primordial non-Gaussianity. Unlike the case with cosmic microwave background, 
21-cm anisotropies provide multi-redshift information with frequency selection and is not damped at 
arcminute angular scales. We discuss the angular trispectrum of the 21-cm background anisotropics 
and discuss how the trispectrum signal generated by the primordial non-Gaussianity can be measured 
with the three-to-one correlator and the corresponding angular power spectrum. We also discuss the 
separation of primordial non-Gaussian information in the trispectrum with that generated by the 
subsequent non-linear gravitational evolution of the density field. While with the angular bispectrum 
of 21-cm anisotropies one can limit the second order corrections to the primordial fluctuations below 
/nl ~ 1, using the trispectrum information we suggest that the third order coupling term, /2 or 
<7nl, can be constrained to be arounde 10 with future 21-cm observations over the redshift interval 
of 50 to 100. 

PACS numbers: 98.80.Es,95.85.Nv,98.35.Ce,98.70.Vc 

I. INTRODUCTION 

The cosmic 21-cm background involving spin- flip line emission or absorption of neutral hydrogen contains unique 
signatures on how the neutral gas evolved from last scattering at z ~ 1100 to complete reionization at z < 10 
Subsequent to recombination, the temperature of neutral gas is coupled to that of the cosmic microwave background 
(CMB). At redshifts below ~ 200 the gas cools adiabatically, its temperature drops below that of the CMB, and 
neutral hydrogen resonantly absorbs CMB flux through the spin- flip transition 0, H[ 0] . The inhomogeneous neutral 
hydrogen density distribution generates anisotropies in the brightness temperature measured relative to the blackbody 
CMB Q. The large cosmological and astrophysical information content in 21-cm background is well understood in 
the literature [1 @ 0, S H M, H3 • 

Parallel to the large effort in analytical and numerical calculations of the 21-cm properties, there are now several 
first generation 21-cm experiments underway focusing on the 21-cm signal during the era of reionization. At the 
low redshifts probed by these first generation interferometers, the 21-cm signal is modified by the astrophysics of 
first sources and the associated UV photon background. There, one naturally expects fluctuations to be dominated 
by inhomogeneities in source properties [T2T [l3j . The associated non-Gaussianity in 21-cm anisotropies leads to a 
measurable three-point correlation function or a bispectrum [l4l flBl fl6| . Such a non-Gaussianity is also expected 
to dominate the signature in 21-cm anisotropies generated by the primordial non-Gaussian density field. However, 
the 21-cm background generated by neutral hydrogen at redshifts of 50 to 100 prior to the onset of reionization 
and the appearance of first stars is expected to provide a cleaner probe of the primordial density perturbations in 
the same manner CMB observations are used to study primordial fluctuations [3| . If the primordial fluctuations are 
non-Gaussian, then the 21-cm anisotropies at these high redshifts will naturally contain a signature associated with 
that non-Gaussianity [13, HI] 

The primordial non-Gaussianity in the density field can be studied with the three-point and higher-order correlation 
functions of the 21-cm background. In particular, the second order corrections to the density perturbations generated 
by primordial non-Gaussianity lead to a bispectrum with a dependence on the second-order correction to the curvature 
perturbations, /nl [HI IHj]. With future low frequency data out to z ~ 100, 21-cm background anisotropies could 
potentially limit /nl < 0.1 17]. In comparison, the expected non-Gaussianity under standard inflation is of order 
1 7i s — 1| and with the scalar spectral index n s ~ 0.98 [29| . /nl is expected to be well below unity [25j]. The primordial 
non-Gaussianity parameter at the second order /nl, however, has a correction associated with evolution of second 
and higher-order perturbations after inflation |26j. For standard slow-roll inflation then /nl = — 5/12 (n s — 1) + 5/6 + 
3/10/ (&) where the last term is momentum dependent. In this case /nl is at the level of a few tenths and could 
be as high as 1. The ability for 21-cm anisotropies to probe /nl as low as 0.1 is imp ortant since even a perfect 
CMB experiment limited by cosmic variance alone can only restrict /nl > 3 [HI, [2(| \2l\ while there is no significant 
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improvement when using low redshift large-scale structure [23l . l24j . 

The possibility to make primordial non-Gaussianity measurements with the 21-cm background is important since 
compared to most other probes of inflationary parameters, /nl is one of the few parameters for which we have 
limited number of probes sensitive to the low amplitude non-Gaussianity expected under standard slow-roll inflation. 
While CMB as a probe of non-Gaussianity is well known, when compared to CMB temperature and polarization 
anisotropics, 21-cm background has two distinct advantages: (1) the ability to probe multiple redshifts based on 
frequency selection and (2) the lack of a damping tail in the 21-cm anisotropy spectrum, unlike damping of CMB 
anisotropics at a multipole around 2000. 

While the 21-cm background as a potentially interesting and a useful probe of /nl is now known, different scenarios 
for primordial fluctuations may produce a small /nl but a large third-order correction. In certain alternative models 
of inflation higher order terms may be significant even if the second-order term is small and such scenarios include 
the new ekpyrotic cosmology [27j and, under certain conditions, the curvaton model [28| . Thus, beyond the non- 
Gaussianity at the three-point level with the bispectrum, it is also useful to study the non-Gaussianities of the 21-cm 
background at the four point-level involving the trispectrum. 

Here we show that the angular trispectrum of 21-cm anisotropies contains a measurable non-Gaussianity from 
primordial fluctuations if the scale-independent cubic corrections to the gravitational potential captured with an am- 
plitude parameter (or <7nl as described in Ref. (30j ) has a value of order ~ 10, even if the scale-independent quadratic 
corrections to the gravitational potential captured by /nl has a value around ~ 1. While /nl is currently constrained 
with WMAP data [29J, there is no real constraint on this third order non-Gaussianity parameter. There are, however, 
theoretical expectations: for slow-roll inflation, a trispectrum of the form T4 — l/2T^i\P(ki)P(k2)P{k^) + ...] is 
expected for curvature perturbations [3l| with an expectation value of r < r/50 where r is the tensor-to-scalar ratio 
(r < 0.6 is recent CMB data [29j]). In such a model, there is also a direct connection between tnl and /nl such 
that tnl = (6/nl/5) 2 (see the discussion in Ref. [32] for connections between coupling terms of various models for 
primordial trispectrum). While such a relation is generally assumed when constraining tnl [33| . in future it may be 
that one can directly test the above relation between /nl and tnl with data. 

To measure the non-Gaussianity at the four-point level, we introduce the three-to-one correlator statistic, extending 
the two-to-one correlator of Ref. [34[ and applied to WMAP data in Ref. [35| . We optimize the angular power spectrum 
of the 3-1 correlator to detect the primordial trispectrum by appropriately filtering 21-cm anisotropy data to remove 
the non-Gaussian confusion generated by the non-linear evolution of gravitational perturbations. To do this properly, 
one requires prior knowledge on the configuration dependence of the primordial 21-cm trispectrum, but its amplitude 
is a free variable to be determined from the data. 

This paper is organized as following: we first discuss the bispectrum in 21-cm anisotropies associated with primordial 
perturbations resulting from quadratic corrections to the primordial potential. We discuss ways to measure this 
bispectrum in the presence of other non-Gaussian signals and determine the extent to which /nl can be measured 
from 21-cm background data. In the numerical calculations described later, we take a fiducial flat-ACDM cosmological 
model with 6 = 0.0418, tt m = 0.24, h = 0.73, t = 0.092, n* = 0.958, and A(k = 0.05 Mpc -1 ) = 2.3 x 10~ 9 . This 
model is consistent with recent measurements from WMAP [291. 



II. CALCULATIONAL METHOD 



The 21-cm anisotropies are observed as a change in the intensity of the CMB due to line emission or absorption at 
an observed frequency v. 

rp , ~ ^ T S -TCMB (« v n , 

' ' ' = — TT^ — T ^ n ' v ' ^ ' 

where Tg is the spin temperature of the neutral gas, z is the redshift corresponding to the frequency of observation 
(1 + z = Vn/v, with v 2 i = 1420 MHz) and T CMB = 2.73(1 + z)K is the CMB temperature at redshift z. The optical 
depth, t, in the hyperfine transition [2j, when accounted for density and velocity perturbations of a patch in the 
neutral hydrogen distribution [1, 0, 0] , is 

_ 3c 3 /L4i n H f 1 , x 1 + z dv \ 
T lQk B vt l T s H{z)V + dH H(z)dr)> ( ' j 

where Aiq is the spontaneous emission coefficient for the transition (2.85 x 10~ 15 s _1 ), nni is the neutral hydrogen 
density, 8h — ("h - ^h) / is the inhomogeneity in the density, v is the peculiar velocity of the neutral gas distribution, 
r is the comoving radial distance, and H(z) is the expansion rate at a redshift of z. For simplicity, we have dropped 
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the dependences in location n. The fluctuations in the CMB brightness temperature is 

ST b = T b 
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where S(z) describes the coupling between fluctuations in the spin temperature and the neutral density distribution [U, 
5_h is fluctuations in the neutral density distribution, and v(r) is the peculiar velocity. In above T b 26.7mKVl + z 



We describe fluctuations in the 21-cm background in terms of multiple moments defined as ai m {y) = 
f dnY* m (h)5T b (h,v) where 
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with Ji{x) = ji(x) — j'/(x). In equation ((4]), D(r) is the growth function of density perturbations, M(k) = 
—3k 2 T(k)/5il rn HQ maps the primordial potential fluctuations to that of the density field. 

For simplicity in notation, hereafter, we will drop the explicit dependence on the frequency, but it should be 
assumed that 21-cm background measurements can be constructed as a function of frequency, and thus as a function 
of redshift, with the width in frequency space primarily limited by the bandwidth of a radio interferometer captured 
by the window function W u (r). 

Similar to calculations related to the CMB anisotropies, we can introduce a transfer function and rewrite multipole 
moments of 21-cm anisotropies as 



a In 



4Tt(-i) 1 



d 3 k 



$ prim (k). 921cnU (fc)^„(k), 



with the transfer function of the 21-cm background anisotropies given by 
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For reference, the angular power spectrum of 21-cm anisotropies follows by taking (ai m af, ,) — Ci5u>8 mm i as 
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where the power spectrum of Newtonian potential is defined such that ($L(k)$2(k')) = (27r) 3 ^(k — k')P$$(fc). 



A. 21-cm Trispectrum 

Using the multiple moments of the 21-cm background a; m , we can construct the angular average trispectrum, which 
is a rotationally invariant correlation function [401 ] - We define this the usual way such that 

(a hmi a hm2 a hm3 a Umi ) = J](-1) M f _^ J f ^ ^ ^ J T$*(L) , (9) 

where T}^(L) is the angular averaged trispectrum. Here ii,/2>^3i^4 form a quadrilateral with L as the length 
of the diagonal and matrices are the Wigner 3-j symbols. These symbols are non-zero under these conditions: 
\h ~ h\ < L < li + 12, \h — h\ < L < I3 + £4, l\ + I2 + L =even, Z3 + 1± + L =evcn, mi + 7712 — M, and + = —M. 
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The angular averaged trispectrum has two parts involving a connected piece associated with non-Gaussianity and 
an disconnected part that is non-zero even if fluctuations are Gaussian with TtJ^L) = T G \ l \l(L) +T NG j 1 j 2 (X) [H- 
The Gaussian disconnected part is 



T G \±{L) = (-1)^+V(2^ + l)(2/ 3 + l)C h C l3 6 hl2 6 l3h 6 LQ 

+ (2L + l)C h C l2 [(-iy^ +L S hl3 S l2h + S hU S l2la ] . (fO) 

The connected part can be simplified based on permutation symmetries as 



T NG \±(L)=P^(L) + (2L + 1)^ 

L' 

where 



( iVa+is/'l l ? L 1 p hl 3(r r\ , / -i \L+L' f W h L \phli( T i\ 



(11) 



P££(L) = T NG \±(L) + (-l) l ^ +l * +l * +2L T NG \f 3 (L) + (-l) l ^+ L T NG \±(L) + (-1 ) h+h+L T^ (L) . (12) 

In above, matrices are now the Wigner 6j symbol and 'Tj' 1 / 2 (L) is what is described as the reduced trispectrum in the 
literature [13, |10,|iH. 

To generate a non-Gaussianity in the 21-cm background that will lead to a trispectrum, we assume quadratic and 
cubic corrections to the Newtonian curvature such that 

$(x) = $ L (x) + /nl [$i(x) - <$|(x))] + / 2 $|(x) (13) 

when $l(x) is the linear and Gaussian perturbation and /nl and fa are the coupling parameters, which may or may 
not be scale dependent. Note that / 2 is also identified in some publications as <?nl [30( though we follow the notation 
of Ref. [Hj]. Under this description, existing WMAP data limits —30 < /nl < 74 at the la confiden ce l3oT| . while with 
an ideal CMB experiment fundamentally limited by cosmic variance one can constrain |/nl| < 3 [19j] , The angular 
bispectrum of 21-cm fluctuations, especially if measured between z ~ 30 and z ~ 75 prior to the formation of first 
sources and out to angular scales corresponding to multipoles of ~ 10 4 can be used to limit |/nl| < 0.1 [171 ] . This 
is well below the standard expectations for /nl even after accounting for second-order evolution during horizon exit 
and re-entry [26] ]. There are no useful observational limits on / 2 and no estimate on how well / 2 can be established 
with the CMB trispectrum. 

In Fourier space, we can decompose equation (| 13[) as 



with 
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Using these, the reduced connected trispectrum can be written as 
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with two terms involving 
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Using the 21-cm transfer function defined in equation (7) and using equation (JTTJ) , we write 
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We simplify further by expanding the 8d functions, for example 



Mk 3 + k4 - K) = J AL e — (k3+k 4 -K) f (20) 

and combining with the Rayleigh expansion of a plane wave 

e* rk = (4tt) £ i l j t (kr)Y; m (k)Y^ (f ) (21) 



I m 

to simplify the reduced trispectrum of 21-cm anisotropies as 
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hl ^ L = V 4^ J ' (23) 

Substituting for the reduced trispectrum of Newtonian curvature, and similar to the CMB trispectrum [32l [4l|. 
the reduced trispectrum of 21-cm anisotropies can be written with two contributions from second- and third-order 
corrections as 

TttG^liL) =4/nl J r\dr x J r 2 dr 2 F L (ri,r2)ai^ri)0i 2 (ri)ai 3 (r 2 )f3 h (r 2 )h h t 2 Lhi 3 i iL 
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B. 3-1 power spectrum 

Instead of measuring the full trispectrum to extract information on the non-Gaussianity of primordial perturbations 
captured by /nl and f 2 , we consider a compact statistic that measures non-Gaussian information but can be described 
as a higher order 2-point statistic. For this purpose, we introduce the three-to-one correlator analogous to the two- 
to-one correlator of Ref . (34[ and applied to limit /nl from WMAP in Ref . [35[ • This statistic is 

W(h,m) = (T b 3 (n)T fc (m)) (26) 

I 1 m 1 l 2 m2 

where af m = J dnT 3 (n)Y, m *(n), where T 6 (n) is the 21 -cm brightness temperature with the filter described below 
applied in multipole space. We note that this statistic has been recently discussed in the context of separating lensing 
and kinetic Sunyaev-Zel'dovich contributions to arcminute scale CMB anisotropies [42l ]. 

Similar to the filtered two-to-one angular power spectrum [l7l |34| , the three-to-one cubic angular power spectrum 

is 

(aL4m<> = Xl ot Su>5 mm > , (27) 
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where 



x ' 0t = T^wTt z2 ~ 7Kr—r^ T t°t\l\ 2 ( L ) w hi2h\i,L h hhLhi 3 iL , (28) 

V2l + 1 hl 2 l 3 L K 2L + L ) 

where w^i^^i^ is the form of the filter function that applies to the cubic field to optimize the detection of 
any particular form of the underlying trispectrum. In real space, this filter can be simply described as af m — 
J dnr b 3 (n)^i(n)W2(n)W3(n)y^ m *(n) where W\ to W3 are real space filters that are applied to the three maps of 
the same field that are multiplied together. In above Xi note that contributions come from both the unconnected 
Gaussian part of the trispectrum T g and the connected non-Gaussian part of the trispectrum T c . For simplicity, 
we identify separately Xf rim and Xf rav as the contributions resulting from the trispectra of the density field gen- 
erated by primordial density perturbations (equation |24|) and the subsequent non-linear gravitational evolution of 
the density field (see below), respectively. Here, we are primarily interested in detecting the non-Gaussian infor- 
mation contained in the non-Gaussian part of the trispectrum produced by primordial non-Gaussianities or Xf mm . 
With the Gaussian contribution to the trispectrum, the total three-to-one angular power spectrum can be written as 
X tot = x plim + x siav + x Gaussian 



I — /l, -r -r ^1 

Substituting Eq. 10 in in Eq. 19 we derive a simplified expression for xp aussian as 



^Gaussian = V% + 1 ^ ^ ^ + ^ 

47T L — ' 

hL 



Wl^m^LO + (2L + 1) ( q o) i W lhl\LL + Wn lh \i, L } 



(29) 



If the filter function u^^^l is designed such that it is equal to zero if any of (Zijfojfa) is equal to another then 

^Gaussian = q &nd there 

is no contribution to the three-to-one angular power spectrum from the Gaussian term. 
The variance of the three-to-one power spectrum calculated from {XiX[) — (Xi) 2 with 

{XiX[) = — ^ (af m ai m af^ai m ) , (30) 



is 



(21 + l)Nf = (Xf rim ) 2 + (Xf rav ) 2 + (x ; Gaussian ) 2 



C, tot ^ wf 



E hhh\l,L ^tot/^tot^itot u2 h 2 Ci-\\ 
(2L + 1) 1 2 3 llhL hlL ' ^ ' 

and we have ignored the covariance generated by non-Gaussian terms involving three- to eight-point correlations of 
aim- 

To optimize the detection of X pnm , under the assumption that Xf lav < Xf llm first, we find the shape of the filter 
that maximizes the signal-to-noise ratio for its detection is 

(_lVi+ i a+ £ 'TP"™' 1 ' 2 

u>hi 2 i 3 \i,L - -t 7 n n n ' ^ ' 

with the additional constraint that tfi 1 ; 2 i 3 u l £ = if two of (Zi, I2, h, I) are equal. 
With the filter applied, the noise related to Xf Ilm is 



iV, tot 



/ ygrav\2 „ tot y.2 

\Al I 1 °Z \ " l x l 2 l 3 \l,L ^,tot^tot^tot;,2 .2 



1 (21 + 1)2 (2L+1) 11 



1/2 

(33) 



With the assumption that Xf lav w then maximum signal-to- noise ratio for a detection of X ; prim with a noise 
spectrum of 7V) tot is 
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FIG. 1: The power spectrum of 21-cm anisotropies (top solid line) and the angular power spectrum of the projected three-to- 
one correlator described in the paper as a probe of the 21 cm trispectrum. In addition to X ; prlm with /nl = 10 and /2 = 1, we 
also show Xf rav and Ni. 



This is the same signal-to-noise ratio for a detection of the trispectrum generated by primordial non-Gaussianity. 
However, Xf Tav is not necessarily zero and with the same filter applied in the presence of non-negligible non-Gaussianity 
from non-linear density perturbations, there is a residual contribution to Xi that reduces the overall signal-to-noise 
ratio to be below that of the maximal value for a detection of the primordial trispectrum alone. In this sense, in 
the presence of secondary non-Gaussian signal, the filter in equation Q32p is non-optimal and could potentially be 
redesigned to improve the overall signal-to- noise ratio for a detection of Xf llm . While we do not make such an 
attempt here, in estimating the signal-to-noise ratio, we do account for the contamination from Xf IW and the overall 
signal-to-noise ratios we calculated from Xf llra and Ni is below the signal-to-noise ratio given in equation (|34[) . 

This degradation in the signal-to-noise ratio comes from the cross-correlation of trispectra of primordial non- 
Gaussianity (as used in the filter) and the non-Gaussianity generated by gravitational evolution of density pertur- 
bations. To understand this confusion associated with non-linear gravitational evolution, we also calculate Xf rav 
following the derivation of Tj 51, as described in the Appendix. Since Wi^i^^l is defined in terms of the pri- 

mordial trispectrum, Xf' av <x J2 hh i 3 L T?™^ \L) x T^\f{L) while Xf™ oc J2 hhh L (£)| a - Since modes 

of T gravW2^ do nQt align with those of T c P rim ;^ 2 (i) the former sum has cancellations and in general we do expect 

Xf ra,v to be at the same order as X prim or below. 

The dominant contribution to A^ ; tot is not X ; grav but is the Gaussian variance captured by the second term of 
equation (|33[) . This statement is independent of /nl and ji. If stated differently, when properly filtered to search 
for the primordial non-Gaussianity, the main confusion for detecting primordial signal is not the non-Gaussianity 
generated by non-linear perturbations but rather the Gaussian covariance associated with the statistical measurement 
of Xf llm . In practice, the measurement of Xf is likely to be further confused by the non-Gaussianity of foregrounds, 
which we have mostly ignored in the present discussion. Unfortunately, little is known about the expected level of 
the foreground intensity in the low radio frequency range of interest. Techniques have been suggested and discussed 
to remove foregrounds below the detector noise levels (38l . l39j and the filtering process we have outlined will further 
reduce the confusion from the remaining residual foregrounds. This is clearly a topic for further study once data 
become available with first-generation interferometers [471 ] . 
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III. RESULTS AND DISCUSSION 

In Fig. 1 we summarize the power-spectrum of 21-cm anisotropics generated by the neutral hydrogen distribution 
at a redshift of 100 with a bandwidth for observations of 1 MHz. Here, we also plot X prim , Af rav , and Nf ot for 
the same redshift with the optimal filter applied with /nl = 10 and fa = 1 as the non-Gaussian scale-independent 
amplitude of the primordial second- and third-order curvature perturbations, respectively. As shown, Nf ot > A grav , 
suggesting that the noise term is dominated by the Gaussian variance (second term in equation [33]). This statement is 
independent of /nl and fa and thus the non-Gaussian detection is dominated by the Gaussian term in Ni regardless 
of what is assumed about non-Gaussianity. Note that we have estimated Nf ot in Fig. 1 under the assumption that 
observations are limited only by the cosmic variance and not accounting for any instrumental noise variance, which 
will also lead to a cut-off in I out to which we can make measurements. Using the cosmic variance alone allows us 
to establish the potentially achievable limit and compare directly with cosmic variance limit with CMB data. When 
calculating A prim and A grav in equation (f2"5)) we set maximum value of L in the sum associated with T lll2l -s' (L) to be 
L max = 100 We tested our calculation for L max = 150 and found results to be within a percent, but such a higher value 
slows the numerical calculation significantly. Finally, due to computational limitations of the numerical calculation, 
we restrict estimate of X[ to I = 10 3 . 

In Fig. 1 when calculating A glav , following the derivation in the Appendix and the discussion there, we use the 
exact analytical result for the 2-2 trispectrum of the density field with the mode coupling captured by F^ki, k^) 
[43|,[44[. For the 3-1 trispectrum of the density field under non-linear gravitational evolution, given that an analytical 
result for the angular trispectrum is cumbersome, we use the angular averaged value for F^ki, k2, ^3). We refer the 
reader to the Appendix for details. 

In Fig. 2 we summarize the estimate related to signal-to-noise ratio for a detection of X pnm as a function of /. 
The typical signal-to-noise ratio, when measurements are out to a multipole of 10 3 , is at the level of ~ 0.5 if one 
assumes that the coupling parameters /nl = 10 and fa = 3 for 21-cm observations centered at a redshift of 100 over a 
bandwidth of 1 MHz. The signal-to-noise ratio for the case with /nl = and fa — 3 is ~ 0.03. Since S/N oc fa when 
/nl = 0, out to / = 10 3 , a signal-to-noise ratio of 1 is achieved if fa ~ 10 2 . While there are neither strong theoretical 
motivations on the expected value of fa nor a real bound on its value from existing data, it is likely that with 21-cm 
data one can constrain fa down to a level well below this value for a single redshift. This is due to the fact that 
21-cm observations lead to measurements at multiple redshifts, though one cannot make arbitrarily small bandwidths 
to improve the detection since at scales below a few Mpc, anisotropies in one redshift bin will be correlated with 
those in adjacent bins [3{J. For example, if 21-cm observations are separated to 30 independent bins over the redshift 
interval 50 to 100 (as can be achieved with 1 MHz bandwidths), then an approximate estimate of the cumulative 
signal-to-noise ratio, S/N = ^/J^[S/N(z)] 2 , is ~ 0.15(/2/3) if /nl = 0. In return, one can potentially probe /nl 
values as low as 20 roughly out to Z max ~ 10 3 . With the first-generation radio interferometers, we would at most 
survey 1% of the sky. Assuming instrument noise is dominating at multipoles above 10 3 between 30 MHz and 60 
MHz at 1 MHz bandwidths (corresponding to redshifts 30 to 100), we find a signal-to-noise ratio of 5 x 10~ 3 /2, which 
could lead to a limit on fa of order 200. 

Above discussion on the application of X? nm assumes that /nl = 0. Since X prim oc / NL , if /nl is greater than 
one, the overall signal-to-noise ratio for the detection of the three-to-one angular power spectrum is increased and the 
dominant contribution to A prim comes from the coupling associated with /nl and not fa. In the case when fa = 0, 
the /nl one probes with the trispectrum can be related to tnl of Refs. [31], [33| for the primordial trispectrum. In 
Fig. 2, we show the signal-to-noise ratio with /nl = 10 and fa = 3. Since in this case /nl term dominates, this 
provides an approximate estimate of the signal-to-noise ratio for /nl with the trispectrum. Using a single redshift 
bin out to ^max ~ 10 3 , 21-cm observations achieve a signal-to-noise ratio of 1 if /nl ~ 15. This in return constraints 
iaitNL < 300. Assuming 30 redshift bins over the redshift interval 50 to 100, assuming fa — 0, we find that one can 
constrain tnl < 50. This result is only out to i max = 10 3 , but since 21-cm observations are not damped as in the 
case of CMB observations, higher resolution data can improve limits on both fa and tnl significantly especially if 
observations can be pushed to I > 10 4 . 

While a detection of the CMB angular trispectrum has been motivated as a way to constrain /nl or tnl [3^ |. 
this is probably not necessary with 21-cm data. Once radio interferometers start probing the redshift interval of 50 
to 100, the angular bispectrum, which can be probed with the two-to-one angular power spectrum (l7j . can limit 
/nl < 0.1. This will facilitate a separation of the contribution to the three-to-one power spectrum from /nl and fa- 
While we have assumed that /nl and fa are momentum independent, it is likely that for specific models of inflation, 
these coupling terms are momentum dependent [46| and then the ability to separate /nl and fa by combining 21-cm 
bispectrum and trispectrum information with the two-to-one and three-to-one correlator respectively will strongly 
depend on the exact momentum dependence of the coupling factors. We leave such a study for future research. While 
measuring fa is well motivated, one can also test the consistency between /nl probed by the bispectrum and the 
tnl from the trispectrum related to the slow-roll inflation predictions for the non-Gaussianity. While we have not 
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10 1 p ^ ^ 




FIG. 2: Signal-to-noise ratio for a detection of Xf rlra with /nl = 10 and /a = 3 (solid line) and /nl = and /2 = 3 (dot-dashed 
line). For reference, we show the signal-to-noise ratio associated with the detection of the full 21-cm bispectrum with a dashed 
line. 



discussed the extent to which this consistency relation can be established with upcoming experiments after taking 
into account of instrumental noise it will also be useful to return to such a calculation in the future once 21-cm 
observations begin to probe the universe at z > 10. 



IV. SUMMARY 



The 21-cm anisotropics from the neutral hydrogen distribution prior to the era of reionization is expected to be more 
sensitive to primordial non-Gaussianity than the cosmic microwave background due to both the three-dimensional 
nature of the 21-cm signal and the lack of a damping tail at arcminute angular scales. Previous calculations have 
discussed the extent to which 21-cm bispectrum can be used as a probe of primordial non-Gaussianity at the two-point 
level with an non-Gaussianity parameter /nl [l7l fl8| . 

Here, we extend these calculations to discuss the possibility to use a four point statistic of the 21-cm background as a 
probe of the primordial non-Gaussianity associated with the trispectrum. We have calculated the angular trispectrum 
of the 21-cm background anisotropics and have introduced th three-to-one correlator and the corresponding angular 
power spectrum as a probe of the primordial trispectrum captured by both /nl and the third order non-Gaussianity 
parameter /2 (described in some publications as <?nl)- Since the primordial non-Gaussianity is confused with a non- 
Gaussian signal in the 21-cm background generated by the non-linear evolution of the density perturbations under 
gravitational evolution, we have discussed way to separate the two using an optimal filter. While with the angular 
bispectrum of 21-cm anisotropics one can limit the second order corrections to the primordial fluctuations as low as 
/nl ~ 0.1 below the value of ~ 1 expected for inflationary models, using the trispectrum information we suggest that 
one can constrain third order coupling term ji to about few tens. If /nl is large, it could potentially be possible to 
test the consistency between /nl from the bispectrum and the slow-roll non-Gaussianity tnl at the four-point level 
with the relation tnl = (6/nl/5) 2 . We hope to return to such a detailed study in the future. 



Acknowledgments 



This work was supported in part by NSF CAREER AST-0645427 at UC Irvine and by the Moore Foundation at 
Caltech. We thank David Lyth for motivating us to calculate the 21 cm trispectrum. 



10 



V. APPENDIX 



Here we discuss the trispectrum from the non-linear density field. The trispcctrum is generated by both second 
and third-order perturbative corrections to the density fluctuations: 

S(k) = S w (k) + S (2) (k) + 6 (3) (k) , (35) 

where <5^^(k) = <5i; n (k) is the linear density perturbation and 

— ^ / -^(27rfS D (k 1 +k 2 - k)^ lin (k 1 )5 lin (k 2 )F 2 (k 1 ,k 2 ) (36) 
(2n) J (27r) 

5&(k) = f £^ f ^(2tt) 3 / ^(27r) 3 fe(k 1 +k 2 +k3-k)5 lin (k 1 )4n(k 2 )^ n (k3)F 3 (k 1 ,k 2 ,k 3 ), 
J (27r) J (2tt) J (27r) 

where F 2 (k x ,k 2 ) and F 3 (ki,k 2 , k 3 ) are derived in Ref. [il-El]. 

The reduced trispectrum of density perturbations (£(ki)(5(k 2 )(5(k3)(5(k4)) can be written in terms of the connected 
piece as 

(5(k 1 )5(k 2 )5(k 3 )5(k 4 )) = (2tt) 3 J d 3 ^fe(k 1 +k 2 +K) ( 5 I? (k3+k4-K)7;-(k 1 ,k 2 ,k3,k4;K) (37) 

with two terms involving 

T 5 (2) (k 1 ,k 2 ,k 3 ,k 4 ;K) = 4F 2 (k 1 ,K)F 2 Ck 3 ,K)P s (K)P s (k 1 )P s (k 3 ) 

r 5 (3) (k 1 ,k 2 ,k 3 ,k 4 ;K) = F 3 (k 2 ,k 3 ,k 4 )P S (k 2 )P S {k 3 )P S (k 4 ) + F 3 (k 2 M,k 4 )Ps(k2)P S (k 1 )P S {k 4 ) , (38) 

where 

2 



5 ki • k 2 ki ■ k 2 2 / ki • k : 



^ k " k2) = 7 + i^f + i<r + Htb i J • (39) 

and ^(ki, k 2 , k 3 ) is derived in the Appendix of Ref. [43j |. 

To calculate the angular trispectrum of 21-cm anisotropies we make use of the 21-cm transfer function defined in 
equation (7) and equation (|37[) to write 

/ \ I a \it -\TU f d3k l f d3k 4 / j3rz 

\ai xmi ai^ m ^ai imi ai imi ) = (Air) \ — \ - — -3 a K 



(2ny J (2nY 

x(27r) 3 5_D(ki + k 2 + K)S D (k 3 + k 4 - K)7^(ki, k 2 , k 3 , k 4 ; K)g 2 i cm ,i 1 (k 1 )g 2 icm,i 2 (k2)g2icm,i 3 (k 3 )g2icmM( k 4) 
xY^faY^faY^faYZ^fa) . (40) 

The angular trispectrum associated with the <5^ 2 ) term can be calculated numerically in a reasonable time, but the 
exact numerical calculation of the angular trispectrum of 21-cm anisotropies associated with 5^ term is slow. Here, 
we outline the analytical derivation of the trispectrum associated with the 5^ term, but for the 5^ term, following an 
approach similar to prior calculations of the non-Gaussianity associated with the trispectrum of gravitational evolution 
p5j . we employ an approximation with the angular averaged value Rb = (F 3 (kx,k 2 ,k 3 )) = 682/189 [44[ and ignore the 
exact mode coupling resulting from the F 3 term. This assumption allows us to write Tg = Rb[P(ki)P(k 2 )P(k 3 ) + ..]. 

Compared to the derivation in Section II where the trispectrum from primordial perturbations involve a coupling 
term which is momentum- independent, the derivation of the trispectrum associated with non-linear evolution with a 
momentum-dependent term is tedious. We take multipolc moments of the F 2 term 

F 2 (k!,k 2 ) = (4tt) F 2A(h,k2)Y Uma {k 1 )Y l * ama (k 2 ) (41) 

and since F 2 involves terms (ki • k 2 ) n from n — 0, 1,2, l a takes the values of 0,1,2. We outline the contribution for 
a a specific combination of (l a , lb) involving the expansion of the two F 2 terms, the trispectrum calculation generally 
involves a term of the form 

(a limi a l2m2 a l3m3 a hmi ) = 4(4.) >] J d 3 K J — 3 J — z (42) 

x J2 Y Um a (k )Y l * ama (K)Y hmb (k 3 ) Y i* bmb (K)F 2) i (fci, K)F 2t i b (k 3 , K)P,j 1 (K)P^(k 1 )P l j 1 (k 3 ) 
^■^+^+K)e>^+^^g 2lcm ^^ 
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which can be simplified with Rayleigh expansion of the plane waves followed by angular integrals to arrive after some 
tedious but straightforward algebra to 

\Oil- l mi a l 2 m 2 0'l 3 m 3 0'l 4 m 4 ) = 

(43) 

4(4tt) 12 (-0' 1+ ' 3 y kfh...J k 2 dk 2 J K 2 dK J rfdn J rldr 2 F 2 ^(k 1 ,K)F 2h (k^K)P^K)P^k 1 )P 4> {k 3 ) 
x ]T (-l) mi+ " i2+m3+Ml+i4 g 2 icmj 1 (feOffaicm.ia (fc 2 )<?2ic m ,i 3 (h)g 21cmA4 (fc 4 ) 

rn a m b L i Mi L 2 M 2 L 3 M 3 L 4 
Xji2(fc2r-l)ju(fc4r2)jL 1 (fcir 1 )j L2 (fc3r2)^3(^r 1 )j L4 (i^r 2 )/lL lil ;>L2i3ib^L3i2Li^L4UL2 

L\ h ! a \ / L 2 l 3 It \ / L 3 l 2 Li \ I L4 Z 4 L 2 
Mi -mi m a I \ M 2 -m 3 m b J \ M 3 -m 2 -Mi J \ M 4 m 4 M 2 



x | dK^; ma (K)^: m6 (K)r i * iM3 (K)r i4 M 4 (K) 

Including the exact form of the mode-coupling, we obtain the reduced angular trispectrum of 21-cm anisotropics 
due to second-order gravitational perturbation evolution as 

T grav ^ (jL ) =4 ^ J kjdh... J k\dk 2 J K 2 dK Jrjdn J r 2 dr 2 P4K)P^(ki)P4k 3 ) £ 



L1L2L3L4 

y-g2i cm A 1 {ki)g 2 i cm ,\ 2 {k 2 )g 2 i cm ^{kz)g 2 i cm ^(k^ 



X 



^1 ^2 i 1 J ^3 ^4 



where 

S(l a Jb) = I £^ ^ I I ^ £" 2 I hL 1 l 1 l a hL 2 l 3 l b hL 3 l 2 L 1 hL 4 l 4 L 2 h>L i l b LhL 3 l a L , (45) 

As mentioned above, the angular trispectrum with term involves a calculation that is numerically slow given 
the mode coupling resulting from a term involving F3 oc l/(ki + k 2 + ks) 2 (48j . Since we are considering the 3-1 angular 
power spectrum, we can ignore the exact momentum dependence of the non-linear gravity and take the angular 
averaged value of The resulting trispectrum in this case takes a simple form similar to that of the primordial 
trispectrum with a momentum-independent coupling term, and we overestimate the covariance between primordial 
and non-linear gravity trispectra to the 3-1 angular power spectrum estimator. We do not reproduce the derivation 
of the angular trispectrum with i?f, given that it is similar to equation (|23[) . 
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